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We show that the self-organized criticality of the Bean's state in each of the grains of a granular 
superconductor results in magneto-thermal oscillations preceding a series of subsequent flux jumps. 
We find that the frequency of these oscillations is proportional to the external magnetic field sweep 
rate H e and is inversely proportional to the square root of the heat capacity. We demonstrate 
experimentally and theoretically the universality of this dependence that is mainly influenced by 
the granularity of the superconductor. 
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The theory of the self-organized criticalityl'i explains 
the dynamics of the nonequilibrium systems similar to 
sandpiles. The magnetic flux dynamics and, in particu- 
lar, the magnetization relaxation process in a supercon- 
ductor with a strong pinning potential for vortices is an 
interesting and challenging field of application for this 
theoryoO. In the scope of the ideas of the self-organized 
criticality the pinned vortices space distribution arises as 
a result of a subsequent local vortex avalanches, i.e., a 
series of small local flux jumps is establishing the crit- 
ical state. The Bean critical state modelQ successfully 
describes the irreversible magnetization in type-II super- 
conductors by introducing the critical current density j c . 
In the framework of the Bean model the value of the slope 
of the stationary magnetic field profile is less or equal to 
fJ'ojc- It makes the spatial distribution of vortices in a 
superconductor with a strong pinning potential similar 
to the sand particles spatial distribution in a sandpildl. 

The stationary critical state becomes unstable under 
certain conditions when the local flux jumps result in r 
global flux jump driving the system to the normal state 
This instability can.-be preceded by a series of magneto- 
thermal oscillaiipasa. These oscillations have been re- 
ported earlierE3i!3 but never have been studied. 

Each of the local vortices avalanches establishing the 
critical state produces a heat pulse and a temperature rise 
in the superconductor. This temperature rise decreases 
the critical current density for a certain time interval and, 
thus, changes the initial conditions for the subsequent 
vortices avalanches. In other words, the heat pulses pro- 
duced by the local vortices avalanches result in a correla- 
tion mechanism specific for the self-organized criticality 
of magnetic flux motion in superconductors with a strong 
pinning potential. In this Letter we demonstrate that 



this mechanism results in the magneto-thermal oscilla- 
tions arising close to the threshold of the superconduct- 
ing state stability. We focus our study on the dependence 
of the frequency of these oscillations on the temperature 
and the magnetic field sweep rate in case of a granular 
superconductor. We point out that the granularity and 
the randomness of the properties of the superconducting 
grains result in universality of this dependence. 

We begin with the theoretical consideration of the crit- 
ical state stability and magneto-thermal oscillations. We 
propose a one-dimensional model of a granular super- 
conductor treating it as a stack of superconducting slabs 
having the width 2bi (i — 1, 2, 3, . . . , N) randomly dis- 
tributed with a certain mean value b. We assume that 
there is no electrical contact between the slabs and there 
is an ideal thermal contact between them. We suppose 
that the external magnetic field H e (t) is parallel to the 
sample surface (H e is parallel to the z-axis) and the mag- 
netic field sweep rate H e is constant. We assume that the 
critical state arises simultaneously in the entire supercon- 
ductor and therefore in each of the slabs the background 
magnetic, Bi(x,t), and electric, Ei(x,t), fields are deter- 
mined by Maxwell equations 
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We suppose also that close to the instability threshold 
most of the slabs are saturated, i.e., the external mag- 
netic field H e is higher than the Bean field B p = fJ®j c b. 

Magneto-thermal oscillations in the critical state arise 
as coupled oscillations of small perturbations of the tem- 
perature, 9, and the electric field, e. The heat diffusion 
and Maxwell equations determine the spatial and tem- 
poral variations of 6 and e, namely, 
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where C is the heat capacity, k is the heat conductivity, 
and 
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In the critical state the current density, j, is, close to 
j c . In this region the j-E curve takes the forir£3 



J=Jc+n ln(|-) 



(5) 



where To is the voltage criterion at which j c is defined, 
ji determines the slope of the j-E curve, and ji <C j c - 

The relation given by Eq. (|5|) was ficst derived in the 
framework of the Anderson-Kim modelijiij considering 
the thermally activated uncorrelated hopping of bundles 
of vortices T. It follows from this model that ji oc T. In 
the framework of the self-organized criticality the value 
of j\ seemSpto be temperature independent in the range 
of T <C T c bD. We will consider the later case as it is. in a 
good agreement with numerous experimental datatU. To 
be more precise, we suppose that the ratio n — j c /ji 3> 1 
is temperature independent at T <C T c . 

Using Eqs. (Q) - (||) we obtain the equations to deter- 
mine 9(x,t) oc exp(7i) and e(x, t) oc cxp(7i) in the form 
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Let us clarify the following calculation qualitatively. 
Suppose that the initial temperature of the superconduc- 
tor To increases by a small perturbation 9q arising due 
to a heat pulse with the energy 5Qq. This temperature 
increase leads to a decrease of the superconducting cur- 
rents. The reduction of these screening currents results 
in an additional flux penetration inside the superconduc- 
tor. This flux motion induces an electric field pertur- 
bation eo producing an additional heat release &Q\, an 
additional temperature rise 0\, and, consequently, an ad- 
ditional reduction of the superconducting currents. At 
certain conditions this process results in an avalanche- 
type increase of the temperature and magnetic flux in 
the superconductor, i.e., in a critical state instability. 

The critical state is stable if the heat release, 5Q, aris- 
ing in the process of electric field and temperature per- 
turbations development is less than the maximum heat 
flux to the coolant. The value of SQ depends on both 
H e and H e for the unsaturated grains and only on H e 
for the saturated grains. In our experiments most of the 
grains are saturated in the magnetic field region corre- 
sponding to the magneto-thermal oscillations. Therefore, 



the heat release in the unsaturated grains, 5Q U , is small 
compared with the heat release in the saturated grains, 
SQ S . However, the term SQ U is the only magnetic field 
dependent term in the heat balance equation and, thus, 
it determines the value of the global flux jump field Hj. 
In other words, the relatively small heat release in the 
unsaturated grains is tuning the superconducting state 
in a granular superconductor to the instability. 

The heat release 5Q depends on the frequency of the 
magneto-thermal oscillations, w = Im7, in both the sat- 
urated and unsaturated grains. Thus, to calculate the 
value of lo we take into account only the dominating heat 
release arising in the saturated grains. In other words, 
the frequency uj is mainly determined by the response of 
the active media of the saturated grains to the small tem- 
perature and electric field perturbations. The vortices 
avalanches establishing the critical state become strongly 
correlated in the vicinity of the instability threshold. 

Under conditions of our experiments the temperature 
9 is practically uniform over the cross-section of the sam- 
ple. Therefore, to solve Eq. (^|) we consider 9 to be con- 
stant. In a saturated slab the background electric field 
Ei = H e x, where x = corresponds to the middle plain. 
In this case Eq. (g) takes the form 
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with the boundary conditions e(±6j) = 0. Note that the 
characteristic space scale of Eq. (N) is given by 
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The magneto-thermal oscillations exist for low values 
of H e where I ^ 6^. It allows to solve Eq. (^) using the 
WKB approximation that results in 
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We integrate now Eq. ( TO o ver the cross-section of the 
superconductor using Eq. (|10[) for the electric field e and 
we end up with the frequency of the magneto-thermal 
oscillations in the form 
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where To is the mean temperature of the oscillations and 
where we determine the parameter T* as 
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Thus, the frequency of the magneto-thermal oscilla- 
tions in a granular superconductor is proportional to the 
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magnetic field sweep rate H e . The value of the parame- 
ter T* ~ T c /n 2 is a certain constant in the temperature 
range T < T c . It follows then from Eq. that the 
ratio /x = u>y/C(To)/H e is independent on To an d H e , 
i.e., it is a constant characterizing the properties of the 
superconductor and its granular structure. 

We perform an experimental study of the magneto- 
thermal oscillations in a textured YBa2Cu307_5 super- 
conductor growja-from the melt and heat treated after 
the preparationEH. First, we do the measurements using 
the original sample (SO) with the size of 8x9x5.5 mm 3 . 
Next, we cut this sample in two approximately equal 
parts (SI and S2) and we measure the magneto-thermal 
oscillations in SI and S2. In this way we check the inde- 
pendence of the oscillations frequency on the size of the 
sample, i.e., its intrinsic origin influenced mainly by the 
granular structure of the superconductor. 

We perform the magnetic characterization by DC and 
AC measurements in a Quantum Design SQUID magne- 
tometer using the sample S3 with a size of 0.3x0.5x0.8 
mm 3 (the sample S3 is cut from the sample S2). We find 
the onset temperature at zero magnetic field T c ps 88 K 
and the AC susceptibility transition width AT ps 5 K. 

We show in Fig. 1 one quarter of the DC magnetiza- 
tion loops measured for the sample S3 for several tem- 
peratures in the magnetic field range < H e < 5 T. We 
find also that in parts of the sample SO smaller than S3 
the Bean penetration field remains the same as for S3. It 
means that the space scale of the screening current loops 
is smaller than the dimensions of the sub-samples, i.e., 
the sample consists of small superconducting grains. We 
show in the inset in Fig. 1 the temperature variation of 
the critical current extracted from the saturation value 
of the magnetization. As can be seen the critical current 
density decreases linearly in the range 2 < T < 6 K. 

The experimental setup for the magneto-thermal os- 
cillations measurements consists essentially of a sample 
holder with a thin hollow powdered graphite column with 
low thermal conductance on which the YBCO sample is 
maintained. A small size (0.2 mm thick, 2 mm 2 surface) 
carbon thermometer is glued to the sample. The entire 
sample holder is maintained under a very low pressure of 
He gas (2x 10 -6 Torr) in order to reduce the heat leak to 
the coolant. The sample is first cooled down to the start- 
ing temperature T s in zero magnetic field. Then the field 
is established at a controlled rate and the temperature 
variation of the sample is measured. 

We measure the magneto-thermal oscillations at tem- 
peratures in the range 2 < To < 7 K and the field sweep 
rates within the interval 15<ff e <60G/s. We show in 
Fig. 2 a typical sample temperature time dependence for 
To = 5 K and H e = 20.4 G/s. At low values of the mag- 
netic field H e (t) there is a slow temperature increase due 
to the small vortices avalanches establishing the critical 
state. Above a certain magnetic field value, the sample 
temperature oscillations appear with a period r in the 
range 10 < r < 70 s and an amplitude increasing in 



time. These magneto-thermal oscillations we mainly ex- 
plore and analyze. A flux jump occurs close to the Bean 
field accompanied by a temperature rise up to about 12 
K with a characteristic time of the order of 1 s. Then the 
sample temperature relaxes to the coolant temperature 
with a rate depending on the heat leak. 

We show in Fig. 3 the magneto-thermal oscillations 
frequency dependence on the field sweep rate, for several 
temperatures and different samples. As can be seen the 
value of lu increases linearly with H e within the accuracy 
of our measurements. 

We show in Fig. 4 the ratio /i = ujy/C(T )/H e as a 
function of the magnetic field sweep rate H e for the sam- 
ples SO, SI, and S2 for different temperatures Ta from 
the interval 2 < T < 7 K. We use the relationliao 

C = 100T + 5T 3 (JKrV- 3 ) (13) 

to calculate the heat capacity and we normalize the val- 
ues of n by the mean value Jl of their distribution. We 
see in Fig. 4 that the ratio fi is a constant within the 
accuracy of our experiments as predicted by Eq. (|Tl|). 
Assuming that T~ x ps n 2 {dM/dT)/M we estimate the 
value of n as n ps 14.5 which ijS-jn a good agreement with 
the known experimental datatij. 

To verify the uniformity of the sample temperature 
we estimate the characteristic heat redistribution time 
th = Cd 2 / k, where d is the sample size. We find that 
t h ps 0.06 s using the data T ps 5 K, C ps 1200 JKTW 3 , 
k ps 2 WK _1 m _1 , and d ps 1 cm. This time constant 
is two orders of magnitude less than both the period of 
the magneto-thermal oscillations and the temperature re- 
laxation time. Thus, the temperature of the sample is 
uniform as suggested above. 

In conclusion, we show theoretically that the self- 
organized criticality of the Bean's state in each of the 
grains of a granular superconductor results in magneto- 
thermal oscillations preceding flux jumping. We study 
these oscillations experimentally in a granular YBCO 
samples at temperatures 2 < T < 7 K and field sweep 
rates 10 < H e < 60 G/s. We find both experimen- 
tally and theoretically that the frequency of the magneto- 
thermal oscillations is proportional to the magnetic field 
sweep rate and is inversely proportional to the square 
root of the heat capacity. We demonstrate the universal- 
ity of this dependence measured for different samples. 
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FIG. 1. One quarter of the magnetization loops for the 
sample S3 for several temperatures. The inset shows the tem- 
perature dependence of the critical current density. 

FIG. 2. The temperature of the sample SI as a function of 
the magnetic field H e for H e = 20.4 G/s and T = 5.0 K. The 
inset shows the magneto-thermal oscillations. 

FIG. 3. The dependencies ui{H e ) for the samples: SO at 
T = 6 K (filled circles), SI at T = 3.8 K (filled triangles) 
and at To = 5.3 K (open circles), and sample S2 at T = 5.2 
K (open squares). The error bars are 8%. 

FIG. 4. The ratio fj, = for the samples SO, SI and S2 as a 
function of the sweep rate for different temperatures To from 
the range 2 < T () < 7 K. The values of fj, are normalized by 
the mean value of their distribution. The error bars are 13%. 
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Fig. 1. L. Legrand, I. Rosenman, R.G. Mints, and G. Collin, 

"Self-Organized Criticality Effect on Stablity: Magneto-Thermal 
Oscillations in a Granular YBCO Superconductor". 




Fig. 2. L. Legrand, I. Rosenman, R.G. Mints, and G. Collin, 

"Self-Organized Criticality Effect on Stability: Magneto-Thermal 
Oscillations in a Granular YBCO Superconductor". 




Fig. 3. L. Legrand, I. Rosenman, R.G. Mints, and G. Collin, 

"Self-Organized Criticality Effect on Stablity: Magneto-Thermal 
Oscillations in a Granular YBCO Superconductor". 
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Fig. 4. L. Legrand, I. Rosenman, R.G. Mints, and G. Collin, 

"Self-Organized Criticality Effect on Stablity: Magneto-Thermal 
Oscillations in a Granular YBCO Superconductor". 



